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ABSTRACT 

Recent deep Chandra surveys of the Galactic center region have revealed the existence of a faint, hard 
X-ray source population. While the nature of this population is unknown, it is likely that several types 
of stellar objects contribute. For sources involving binary systems, accreting white dwarfs and accreting 
neutron stars with main sequence companions have been proposed. Among the accreting neutron star 
systems, previous studies have focused on stellar wind-fed sources. In this paper, we point out that 
binary systems in which mass transfer occurs via Roche lobe overflow (RLOF) can also contribute to 
this X-ray source population. 

A binary population synthesis study of the Galactic center region has been carried out, and it is 
found that evolutionary channels for neutron star formation involving the accretion induced collapse of 
a massive ONeMg white dwarf, in addition to the core collapse of massive stars, can contribute to this 
population. The RLOF systems would appear as transients with quiescent luminosities, above 2 keV, in 
the range from 10'^^ — 10^^ ergs s~^. The results reveal that RLOF systems primarily contribute to the 
faint X-ray source population in the Muno et al. (2003) survey and wind-fed systems can contribute to 
the less sensitive Wang et al. (2002) survey. However, our results suggest that accreting neutron star 
systems are not likely to be the major contributor to the faint X-ray source population in the Galactic 
center. 

Subject headings: binaries: close — stars: evolution, formation, neutron — X-rays: binaries 



1. INTRODUCTION 

Recent Chandra observations of the Galactic center 
region (e.g., Wang, Gotthelf & Lang 2002, hereinafter 
WGL02; Muno et al. 2003, hereinafter MM03) have led to 
a renewed interest in the understanding of compact stellar 
X-ray sources. Both of these long exposure surveys have 
revealed a great number of X-ray point sources with lumi- 
nosities, <j 10'^^ ergs s~^ and spectra which can be fit 
by an absorbed power law. 

Although the nature of these low luminosity sources is 
unknown, several possible interpretations have been pro- 
posed. In particular, Muno et al. (2004) suggest that the 
intermediate polar subclass of magnetic cataclysmic vari- 
ables may dominate the faint population characterized by 
luminosities less than lO'^'^ ergs s~^ and a photon power 
law index ranging between and 1 . Among models involv- 
ing accreting neutron stars (NS), Pfahl, Rappaport & Pod- 
siadlowski (2002, hereinafter PRP02) explored the hypoth- 
esis that wind-fed neutron star systems with intermediate 
and massive main sequence (MS) donors (M ^ 3 Mq) may 
account for the majority of sources with = 10^'^ — 10^^ 
ergs s~^. Furthermore, the contributions from a popula- 
tion of pre-low mass X-ray binaries in which the neutron 
stars accrete from the wind of low mass MS companions 
(< 2 Mq) has also been studied by Bleach (2002) and 
Willems & Kolb (2003). However, it was recently pointed 
out by Popov (2003) that the neutron stars in this lower 
mass binary population may not necessarily accrete all 
the transferred material, implying that only a fraction of 
these types of systems would be continuous X-ray emitters. 
Among the non stellar models that have been proposed. 



Bykov (2003) has suggested that the faint hard sources re- 
sult from the interaction of fast moving supernova ejecta 
fragments with the dense interstellar medium of the Galac- 
tic center region. These sources, due to their low intrinsic 
brightness, may only contribute to the deeper MM03 sur- 
vey. 

The aforementioned stellar studies focused on only 
wind-fed NS systems, however there are a number of 
RLOF transient systems which, in quiescence, emit at X- 
ray luminosities very similar to those of the many point 
sources observed in the Galactic center. Specifically, tran- 
sients with neutron stars in their quiescent state are ob- 
served to have X-ray luminosities in the range of 10^^ — 10^'^ 
ergs s""'^ (see Campana & Stella 2003), while transient sys- 
tems with black hole (BH) accretors tend to be observed at 
lower luminosities of ~ 2 x 10'^^ ergs s~^ (e.g., Tomsick et 
al. 2003) in their inactive state. Although many of the NS 
transient systems are characterized by soft spectra in their 
quiescent state and are difficult to detect from the Galactic 
center region, some fraction have a power law component 
extending beyond 2 keV with luminosities in the range of 
~ 10'^^ — 10^^ ergs s~^ (see Jonker, Wijnands, & van der 
Klis 2004). Such systems may have been detected in the 
deeper survey of MM03. The systems with BH accretors 
are also known to emit at energies greater than 2 keV and 
may also contribute to the faint source population at this 
luminosity level. 

X-ray transient systems with large outburst to quies- 
cent flux ratios are known to have either a NS or BH ac- 
cretor and a predominantly low mass companion (classi- 
fied as a low mass X-ray binary, LMXB). The outburst 
phenomenon is generally believed to result from a ther- 
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mal/ viscous instability in an accretion disk surrounding 
the compact object. Only those systems with mass trans- 
fer rates below a critical value (e.g., Dubus et al. 1999; 
Mcnou, Pcrna & Hcrnquist 2002) can exhibit such tran- 
sient behavior. In general, the donors of LMXBs are be- 
lieved to be either low mass MS or giant-like stars. How- 
ever, observations of ultra-short period LMXBs (with or- 
bital periods less than about 80 minutes), so called ul- 
tracompacts, have revealed that, at least, some of them 
must host very low mass (0.01 — 0.1 Mq) white dwarf 
(WD) donors (e.g., Markwardt ct al. 2002, Galloway ct 
al. 2002; see also Deloye & Bildsten 2003). In a recent 
study Belczynski & Taam (2004) showed that the popula- 
tion of LMXBs with low mass WD donors (and thus with 
very low mass transfer rates leading to transient behav- 
ior) may be more significant than previously anticipated. 
Since the majority of transient sources are expected to 
be in quiescence due to the very low X-ray duty cycle, 
wc siiggcst that a fraction of the low luminosity sources 
observed in the Galactic center region may reflect a pop- 
ulation of transient systems in their quiescent state, while 
the same systems in outburst may contribute to a bright 
source population. 

A binary population synthesis model StarTrack (pre- 
sented in §2) is used to determine the number of X- 
ray binaries (NS or BH accretor, wind-fed or RLOF) in 
the Galactic center region. The entire mass spectrum of 
evolved and unevolvcd donors is considered. In particular, 
low, intermediate, and high-mass MS stars are investigated 
to account for the systems considered in previous studies. 
The results of our studies are described in § 3 and their 
implications are discussed in the final section. 

2. MODEL DESCRIPTION 

The StarTrack population synthesis code developed 
originally for the modeling of binaries with two compact 
objects (Belczynski, Kalogera & Bulik 2002) is used. This 
code has recently undergone major revisions (Belczynski 
et al. 2004, in preparation), reflecting the new features 
associated with the improved treatments of the physical 
processes important for the formation and evolution of X- 
ray binary systems. Specifically, the additions and revi- 
sions include a detailed treatment of tidal synchroniza- 
tion and circularization, individual treatment of various 
RLOF mass transfer phases, and a full numerical orbit 
evolution with angular momentum losses due to magnetic 
braking (MB), gravitational radiation (GR), mass trans- 
fer/loss and tides. 

In contrast to population synthesis studies of other 
groups, wc also include the possibility that neutron star 
formation can occur via an accretion induced collapse 
(AIC) of a massive ONeMg white dwarf (Belczynski & 
Taam 2004). Such white dwarfs in binaries form from 
stars within the mass range of ~ 7 — 12 M© (Belczynski 
& Taam 2004; Tauris & van den Heuvel 2003; see also 
Nomoto & Hashimoto 1988) which is wider than for the 
formation of ONeMg cores in single stars (~ 8 — 10 M©) 
due to tidal mass loss (Nomoto & Kondo 1991) and reju- 
vention (Belczynksi & Taam 2004). The ability to accu- 
mulate matter on the white dwarf is enhanced, paradox- 
ically, by the effect of an optically thick wind from the 
white dwarf surface, which can stabilize the mass trans- 
fer in the system at high mass transfer rates (see Kato & 



Hachisu 1994, Hachisu, Kato, & Nomoto 1999). The ac- 
cumulation ratio of hydrogen-rich and helium-rich matter 
is taken from Hachisu, Kato, & Nomoto (1999) and Kato 
& Hachisu (1999) respectively (sec also Ivanova & Taam 
2004). For the direct accretion of helium or carbon/oxygen 
matter onto the ONeMg white dwarfs we make use of the 
work of Kawai, Saio, & Nomoto (1987) in determining the 
evolution of the accreting white dwarf. Our population 
synthesis study of the Galactic center should, therefore, 
be considered complementary to that of PRP02 who did 
not consider the AIC process nor RLOF X-ray sources. 

The tidal effects have been calibrated using the observa- 
tions of binary periods and eccentricities in several stellar 
clusters (Mathicu ct al. 1992), recovering the orbital decay 
rate in high- mass X-ray binaries (Levine, Rappaport & Zo- 
jchcski 2000). The StarTrack RLOF calculation involves 
the detailed calculation of the mass transfer rates based on 
radius-mass exponents calculated both for the donor stars 
and their Roche lobes. The results of our calculations have 
been compared to a set of published RLOF sequences (e.g.. 
Beer & Podsiadlowski 2002) as well as to calculations ob- 
tained with an updated stellar evolution code (Ivanova et 
al. 2003). Our approach to the mass transfer calculations 
allows for the possibility of (i) conservative versus non- 
conservative RLOF episodes (ii) thermally driven RLOF 
versus nuclear /MB / GR losses driven RLOF and (iii) a sep- 
aration of systems as persistent or transient depending on 
whether the donor RLOF mass transfer rate lies below the 
critical rate for instability to develop in the accretion disk 
(adopted from Dubus et al. 1999 or Menou et al. 2002 
for different compositions of transferred material). Per- 
sistent X-ray emitters remain at a X-ray himinosity level 
corresponding to the secular RLOF mass transfer rate. 
Transient X-ray systems remain at very low luminosities 
for most of their lifetime, however, they exhibit outbursts 
with luminosities that can reach the Eddington limit. We 
account for wind accretion onto compact objects follow- 
ing Hurley, Tout & Pols (2002) with the wind mass loss 
rates adopted from Hurley, Pols & Tout (2000) extended 
to include low mass and intermediate mass MS stars uti- 
lizing the formulae of Nieuwenliuijzen & de Jager (1990). 
The wind velocities are assumed to be proportional to the 
escape velocity from the stellar surface with a normaliza- 
tion proposed by Hurley et al. (2002). For all wind- fed 
systems, we assiimc an orbit-averaged accretion rate. The 
flaring behavior associated with Be stars and the perias- 
tron passage enhanced mass accretion onto the NS are not 
taken into account. 

A large sample of single (10^) and binary (10^') stars 
in the Galactic center arc evolved for 10 Gyrs assuming a 
constant star formation rate. We choose parameters from 
the reference model in Belczynski et al. (2002) with a 
few differences. Specifically, the maximum NS mass is set 
equal to 2 M©, and the most recently inferred natal NS 
kick distribution (for NSs formed by the core collapse of 
massive stars) of Arzoumanian, Chernoff & Cordes (2002) 
is incorporated. For cases in which NS formation occurs 
via the accretion induced collapse evolutionary channel, 
no kicks arc applied (see Belczynski & Taam 2004). How- 
ever, wc calculate one model in which we allow for the full 
kicks imparted to NSs formed through AIC, to assess the 
effect of our no AIC kick assumption. In the primordial 
progenitor population, the primary masses are selected in 
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the range 4—100 Mq (with an initial-mass- function slope 
of -2.7; see Kroupa, Tout & Gilmore 1993), and the sec- 
ondary masses are selected in the range 0.08—100 M0 with 
a flat mass ratio (secondary divided by primary mass) dis- 
tribution. In the rapid phases of binary evolution, the 
accretion rate onto the NS and the BH is limited at the 
maximum Eddington limit with the remaining transferred 
material lost with a specific orbital angular momentum 
of the accretor, but allowing for hypcr-critical accretion 
in the common envelope (CE) phases (e.g.. Blondin 1986; 
Chevalier 1989; Brown 1995). In addition to the dynam- 
ically unstable RLOF mass transfer events leading to a 
CE phase, we also allow for evolution into the CE phase 
in the cases in which the trapping radius of the accretion 
flow exceeds the Roche lobe radius of the accretor (e.g., 
King & Begelman 1999; Ivanova et al. 2003). Finally, 
in comparison with our previous study, we have assumed 
that all systems that enter a CE phase with the donor in 
the Hcrtzsprung gap merge into a single star (Ivanova & 
Taam 2004; Belczynski & Taam 2004). 

The updated version of the code has been tested and 
used for a study of Galactic ultra compact binaries (Bel- 
czynski & Taam 2004), Galactic young black hole pop- 
ulations (Belczynski, Sadowski & Rasio 2004) and used 
to reproduce the X-ray luminosity function of the nearby 
starburst galaxy NGC 1569 (Belczynski et al. 2004). 

3. RESULTS 

3.1. Source Type and X-ray Luminosity 

In our study only accreting binaries with NS and BH 
primaries, which are brighter than Lx = 10^° ergs s~^ are 
considered. The secondaries in these X-ray binaries may 
lose material either through a stellar wind or via RLOF. 
In the latter case, the donors transfer all the material to- 
ward the accretor, whereas for the wind-fed systems only 
a fraction of the material is captured by the compact ob- 
ject. For every accreting system the bolometric luminosity 
(Z/boi) is calculated based on the secular average mass ac- 
cretion rate (Mace) as 



where G is gravitational constant. Afacc is the mass of the 
accretor, i?acc is the radius of the accretor (10 km for a NS 
and 3 Schwarzschild radii for a BH), and e gives the con- 
version efiiciency of gravitational binding energy to radia- 
tion associated with accretion onto a NS (surface accretion 
e = 1.0) and onto a BH (disk accretion e = 0.5). 

To determine the simulated X-ray luminosities of our 
population synthesis, we note that very little of the trans- 
ferred material is accreted for transient systems during 
their inactive state. The spectra of neutron stars during 
this state consist of a soft blackbody component and a 
hard power law component with a tendency for a larger 
contribution from the power law component for sources 
characterized by lower quiescent luminosities (Jonkcr et 
al. 2004, see also Tomsick et al. 2004). The soft com- 
ponent is described by energies of ^ 0.1 — 0.3 keV (see 
Campana & Stella 2003). The theoretical models based 
on deep crustal heating associated with pycnonuclear re- 
actions (Brown, Bildsten, & Rutledge 1998; Colpi et al. 



2001) can provide an understanding of the black body-like 
nature of the spectrum and the observed emission level 
(~ 10^^ — 10^^ ergs s~^) without necessarily invoking the 
accretion of matter. However, short term variability dur- 
ing the quiescent state in Aql X-1 (Rutledge et al. 2002) 
and in Gen X-4 (Campana et al. 2004) has also been ob- 
served, suggesting that accretion of matter may be occur- 
ring in Aql X-1 and Gen X-4 (Asai et al. 1996, 1998; Cam- 
pana et al. 1998, 2000; Campana 2004). The accretion 
rate during this phase is likely to be at a low level since the 
inferred himinositics above 2 keV are ~ lO^i erg s ^. In 
contrast to the majority of transient systems in quiescence, 
the spectral energy distribution of the millisecond accret- 
ing X-ray pulsar, SAX J1808. 4-3658, is dominated by a 
hard component with a luminosity of ^ 5 x lO''^ ergs s~^. 
Such a nonthermal component may arise from synchrotron 
emission in a shock front produced by the interaction of a 
relativistic pulsar wind with matter from the companion 
star (Tavani & Arons 1997; Campana et al. 1998). Based 
on these observational results, we adopt 10^^ ergs s~^ as 
a lower limit for the hard X-ray luminosity, above 2 keV, 
for transient NS systems in quiescence, recognizing that 
the average luminosity level can be higher ^ 10"^^ ergs s^^ 
(e.g., Rutledge et al. 2001; Rutledge 2002; see Jonker et 
al. 2004). Due to the lack of a definitive theory for the 
hard X-ray emission component of transient systems in 
quiescent, we take a semi empirical approach for the qui- 
escent NS transients and adopt an X-ray luminosity level 
of 10^^ — 10^^ ergs s~^ above 2 keV. Furthermore, we as- 
sume that the quiescent NS transient X-ray luminosities 
are evenly distributed in the above range. 

The quiescent emission from BH transient systems, on 
the other hand, is likely related to a low level of mass ac- 
cretion. In contrast to the NS transient systems, the spec- 
tra are typically harder and arc not described by a black 
body. Observations of the 14 black hole transient systems 
for which the quiescent emission has been obtained reveal 
luminosities in the range from ~ 10'^° — 10^^ ergs s~^ with 
a median luminosity of ~ 2 x 10^^ ergs s~^ (see Tomsick 
et al. 2003). As for the NS transient systems a semi em- 
pirical approach is adopted for the BH transient systems. 
In particular, we assume that most (80%) of the quiescent 
BH transient X-ray luminosities above 2 keV are evenly 
distributed in the 10^*^ — 10^^ ergs s~^ range, while the 
rest (20%) of the systems are slightly brighter: luminosi- 
ties evenly distributed in the ~ 10^^ — 10^'^ ergs s~^ range 
(see Fig. 3 of Tomsick et al. 2003). 

Transient RLOF systems in outburst are bright X- 
ray sources, with luminosities corresponding to a fraction 
{flout) of^ the critical Eddington luminosity. The long pe- 
riod systems, with orbits that are sufficiently extensive for 
a large accretion disk to be formed, are usually found to 
emit at the Eddington luminosity (icdd) during outburst, 
while the outburst luminosities of short period systems 
are lower by about an order of magnitude. Therefore, we 
apply an X-ray luminosity correction factor at outburst 
corresponding to Tjout =0.1 and rjout = 1 for the short and 
long period systems respectively. The critical periods, over 
which the Eddington luminosity is adopted, are taken to 
be 1 day and 10 hrs for NS and BH transients in outburst 
respectively (Chen, Shrader & Livio 1997; Garcia et al. 
2003; T. Maccarone, private communication). 
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For the persistent sources (both RLOF and wind-fed "'^) 
and all transients in the outburst stage, where accretion is 
the dominant contributor to the observed luminosity, we 
apply a correction factor, T7boi- to the bolometric luminos- 
ity to account for the fraction of energy emitted above 2 
keV. For different types of systems one may expect cor- 
rections of the order of 0.1-0.5. However, for clarity and 
comparison with previous studies the same correction is 
applied to all the sources within the limits T/boi = 0.1 — 1.0, 
with r/boi = 0.1 chosen as a standard value. 

Therefore, we can obtain the simulated X-ray luminosity 
([ergs s"'^]) from: 



1031 _ 1032 
1030 _ 1032 

1032 _ ^o-''^^ 
'?bol'7out-^^Edd 
'Tbol-t'bol 



all quiescent NS transients 
80% quiescent BH transients 
20% quiescent BH transients 
outburst NS/BH transients 
persistent (RLOF and Wind fed) 

(2) 

where LEdd represents the Eddington luminosity. 



3.2. Standard Model Calculations 

The simulated X-ray luminosity distributions of both 
wind- fed and RLOF- fed systems with NS and BH accrctors 
are illustrated in Figure 1. The RLOF- fed systems dom- 
inate the overall X-ray population at luminosities greater 
than Lx = 10^° ergs s~^ with most of the systems lying 
within the range of ~ 10^^ — 10^^ ergs s~^. The major- 
ity of these RLOF systems arc transient X-ray sources in 
quiescence formed through the AIC scenario presented in 
Belczynski & Taam (2004). The number of X-ray systems 
(435 of RLOF-fed and 56 of wind-fed sources) were cali- 
brated to represent the Galactic center population in the 
MM03 survey and are also listed in Table 1 (for details see 
§3.4). 

In Table 1 the relative frequencies for various types of 

systems categorized by their donor and accretor found in 
our simulation are listed. As expected from the examina- 
tion of Figure 1, the RLOF-fed systems dominate (88.6%) 
the population. The majority of these systems are NS 
transient sources in quiescence with typical luminosities of 
- 10^1 - 10^2 ergs s"!. Most of tlic RLOF-fed systems 
have low-mass WD donors with the mass transfer rates of 
~ 10^^^ M0 yr^^ leading to transient behavior through a 
long lived (~ 1 Gyr) phase (for details see Belczynski & 
Taam 2004). Over half of the systems belong to the ultra- 
short period LMXB class with orbital periods below 2 hrs 
(60.7%), while the remaining RLOF-fed binaries (27.9%) 
are characterized by longer periods extending beyond one 
month. 

A small, but significant fraction of wind-fed systems 
(11.4%) are found, mostly with low-mass (M < 3 Mq: 
7.9%) and intermediate-mass companions (3 < M < 
8 Mq-. 3.3%). However, there is no statistically significant 
population of systems found with massive (M > 8 M0) 
companions (HMXBs). Most of the wind- fed systems are 
NS-MS binaries, liowc;vcr several have an evolved (red gi- 
ant) donor, and a few have a BH accretor. For the wind-fed 
systems, most are characterized by eccentric orbits with 
the most massive donors possibly exhibiting a Be X-ray 
binary phase, although we do not follow orbit-dependent 
accretion. 



Our X-ray population will remain confined to the Galac- 
tic center region since the spatial velocities induced as a 
result of the mass loss and nascent kick velocities accom- 
panying NS/BH formation arc small compared to their 
mean spatial velocities in this region. In particular, the 
average velocities acquired are ~ 60—100 km s~^ for both 
the wind-fed and RLOF mass transferring systems. They 
are significantly smaller than the velocities of ~ 200 km 
characteristic of the stellar motions in this region. We 
point out that the velocities calculated for the RLOF mass 
transferring systems only refer to a minor subgroup (71 out 
of 435 systems, see Table 1) that have undergone a core 
collapse supernova since the majority of these systems are 
formed via the AIC mechanism with no nascent kick. The 
low velocities are a direct consequence of the use of the 
recent NS kick distribution of Arzoumanian et al. (2002). 

3.3. Parameter study 

To determine the sensitivity of the results to the major 
uncertainties of the binary evolution, we have recalculated 
several models with different choices of evolutionary pa- 
rameters. In particular, a model with a reduced efficiency 
of the CE ejection of a x A = 0.1 (in the standard model 
a X X — 1) was simulated; for details see Belczynski et al. 
(2002). In this case, the number of NS/BH X-ray binaries 
(XRBs) in the Galactic center is significantly decreased 
(109), with the relative number of wind- fed to RLOF- 
fed systems remaining comparable to the standard model. 
The sensitivity of the model simulation to the assumed 
mass ratio distribution was also studied. Specifically, a 
model in which both masses were drawn independently 
(as opposed to c:orrclatcd initial masses for the two binary 
components via a fiat mass ratio distribution in the stan- 
dard model) was calculated. This led to an increase in the 
number of systems which enter and merge in the CE phase, 
decreasing the number of NS/BH XRBs in the Galactic 
center (296) and lowering the importance of RLOF-fed 
systems relative to wind-fccl systems (245 RLOF-fed sys- 
tems compared to 51 wind- fed systems). This decrease in 
the number of XRBs for the above two models provides an 
indication of the influence of the uncertainties associated 
with our treatment of the CE evolution phase (for details 
see Belczynski & Taam 2004). 

Since the stellar wind mass loss rates are somewhat im- 
certain, an additional model with the stellar wind mass 
loss rates reduced by a factor of 2 for all mass losing stars 
was calculated. The smaller wind mass loss rates have two 
opposing tendencies on wind- fed systems. First, there is 
an obvious decrease of the mass accretion rate, and thus a 
decrease in the number of wind-fed X-ray systems above 
certain critical value of L^- Second, the systems tend to 
be tighter, due to the fact that less mass is carried away 
(Jeans mode of mass loss). Therefore, the fraction of mat- 
ter lost in the wind that is captured increass, leading to a 
tendency for an increase in the accretion rate. The latter 
effect dominates over the former, and there is an overall in- 
crease in the number of wind-fed systems (82) as compared 
to the standard model (56). The number of RLOF-fed sys- 
tems is not greatly affected (416). 

A model in which we use the full natal kicks (Arzou- 
manian et al. 2002) for NSs formed via AIC as opposed 



^Note that all of the wind- fed systems are modeled here as persistent X-ray sources. 
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to no AIC kicks in the standard model was also calcu- 
lated. In this model wind-fed systems are not affected, 
however a number of RLOF systems progenitors are dis- 
rupted at the formation of the NS. Additionally, many 
of the surviving RLOF systems formed in AIC acquire 
high (> 200km s~^) systemic speeds, exceeding the escape 
velocity from the Galactic center (but still bound to the 
Galaxy). This follows from the result that the AIC chan- 
nel occurs, in general, for systems at short orbital periods 
(where the components orbit about their common center 
at high velocities), and the systems which survive can ac- 
quire a significant additional systemic speed (comparable 
to the component's orbital velocity). The total number of 
NS/BH XRBs decreases to 201, while the relative contri- 
bution of wind-fed sources in the X-ray population signif- 
icantly increases (see Fig. 3) as compared to the standard 
model. 

An additional model, without tidal binary component 
interactions, was calculated to compare the results to the 
other studies. In this model, we find the highest number 
of wind-fed sources (87) in the simulations and a slight de- 
crease of RLOF-fcd systems (387) in comparison with our 
standard model calculation. In the absence of tidal inter- 
action, the orbital expansion associated with transfer of 
spin angular momentum to the orbit does not take place 
and the X-ray binary progenitors are generally found in 
shorter period orbits. For the tightest systems (progeni- 
tors of RLOF sources) this leads to more frequent mergers 
in CE episodes (decreasing the formation rates). On the 
other hand, the absence of tidal effects leads to tighter or- 
bits for wider systems (the progenitors of wind- fed sources) 
and, therefore, to brighter systems resulting from higher 
wind capture rates, increasing their number in the studied 
sample. In particular, the X-ray luminosities of wind-fed 
sources can exceed 10^^ ergs s~^ for intermediate 

mass X-ray binaries. In this calculation, similar to the 
standard model, the number of HMXBs is negligible, as 
expected for the small sampled volume corresponding to 
the MM03 survey. 

We also relaxed our assumption on the bolometric cor- 
rection, and recalculated the NS/BH system X-ray lumi- 
nosities with ?7boi = LO (no bolometric correction). As 
expected (see eq. 2), only the brightest sources in RLOF 
class (either persistent or transients in outburst) are af- 
fected, with Lx extending to ~ 10^^ ergs s"-'^. Since the 
wind-fed systems are brighter (by the factor of ten), the 
number of wind-fed sources increases to 80 as additional 
systems are brighter than our adopted critical threshold 
luminosity {L^ = 10^° ergs s~^). 

Throughout our study we have allowed for the possibil- 
ity of accretion induced collapse of a heavy WD to a NS. 
The specific arguments and new input physics supporting 
AIC was presented and discussed in detail by Belczynski 
& Taam (2004). As an indication of its importance, we 
examine the implications of its neglect on the results of 
our study. The primary effect is a decrease in the number 
of NS/BH X-ray sources to 126 (as compared with 491 
in standard model, see also Table 1). This decrease re- 
sults in a significant depletion of RLOF population (71), 
while wind-fed sources would remain virtually unchanged 
(55). In this model the majority of RLOF systems are 
found with long periods (45), but some systems are still 
found within ultrashort period class (26). As a result, the 



number of transients and the number of bright sources de- 
crease. The shape of the X-ray luminosity distribution is 
not significantly altered, although the relative contribution 
of wind-fed to RLOF-fed sources is increased. 

Notwithstanding the models with rather extreme as- 
sumptions on accretion induced collapse (full natal kicks 
applied during NS formation in the AIC process and a 
model in which AIC is not allowed) the RLOF-fed NS/BH 
XRBs are found to be more numerous than their wind-fed 
counterparts by more than a factor of 2. For comparison, 
the X-ray luminosity distributions for the most extreme 
model assumptions are shown in Figures 2 and 3. Note 
the change of vertical scale between different plots. For all 
of the models most systems are found at low luminosities 
{Lx ~ 10"^° — 10'^^ ergs s^^), with several bright systems 
usually found at (L^; - 10^*^ - 10^® ergs s^^). Although 
the shape of the distribution is very similar, the number of 
sources and the ratio of wind-fed sources to RLOF sources 
may change significantly from model to model. 

3.4. Expected numbers 

The total mass contained in stars corresponding to our 
simulation with 10^ massive binaries (primary component 
more massive than 4Mq) can be estimated based on the 
three component broken law of Kroupa et al. (1999). By 
extending our simulation mass range to O.O8M0 and in- 
cluding a 50% contribution due to single stars, we estimate 
the total mass to be ^ 1.5 x IO^Mq. As based on the recent 
mass model by Launhardt, Zylka, & Mezger (2002), this 
is comparable to the stellar mass in a cylindrical radius 
of 20 pc and a depth of 440 pc centered on the Galactic 
center (Muno et al. 2004). We note, however, that the 
stellar mass would be overestimated by a factor of 5 if the 
MM03 survey corresponded to a sphere of radius 20 pc. 
Our results (Table 1. also previous subsections) can also 
be applied to the WGL02 survey by taking into account 
that the larger field is about 4 times more massive than 
modeled in our population synthesis. In our standard sim- 
ulation it is found that there are ~ 500 and ~ 2000 NS/BH 
X-ray binaries with luminosities Lx > 10^" ergs (see 
Table 1) for the MM03 and WGL02 fields, respectively. 
However, the majority of the synthetic NS/BH accretion 
sources are expected to be characterized by low X-ray lu- 
minosities {Lx ~ 10^^ — 10^^ ergs s^^) in the hard X-ray 
band (see Figure 1, 2 and 3). We note that many of our 
simulated sources can have a significant soft component 
with higher luminosities (see discussion in the §3.1). 

The majority (~ 500) of the Galactic center sources 
discovered in the WGL02 survey were characterized by 
Lx ^ 10'^'^ — lO'^'"' ergs s^^ and hard spectra. In our stan- 
dard model calculation, there are no statistically signif- 
icant sources in that luminosity range. However, in the 
model with no bolometric correction applied (r]boi = 1-0), 
44 sources in that range are found, all of which are wind- 
fed systems. This simply refiects the fact, that all of the 
wind-fed systems found at L^ ~ 10^^ — 10^^ ergs in 
the standard model, are shifted to the higher luminosities 
with ?7boi = 1.0. In addition, for the models with altered 
assumptions on AIC (either full AIC kicks or AIC not al- 
lowed) the number of wind-fed systems (amounting to the 
same as for the standard model) dominate over RLOF-fed 
binaries in this specific luminosity range. Although these 
systems are promising candidates for the unidentified X- 
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ray population, it is likely that they represent a small frac- 
tion 20%) of the observed sources in WGL02 survey. 
Our results suggest that other sources, and not the accre- 
tion powered NS/BH binaries constitute the majority of 
point sources in the WGL02 survey. 

In the deeper survey of MM03 a faint source population 
of 2079 point sources, characterized by ~ 10^° - 10^^ 
ergs s~^ was discovered. The spectra of these sources 
were fit to an absorbed power law with photon index, F. 
The majority of these sources is described by hard spec- 
tra r < 1 (1427), with an important subpopulation char- 
acterized by softer spectra F > 1 (652^). We note that 
the majority of our synthetic XRBs have luminosities in 
a range very similar to these found for point sources in 
MM03 with a significant number of these systems corre- 
sponding to ultrashort period RLOF-fed NS-WD binaries. 
Due to the very low mass transfer rates in these binaries, 
the sources are expected to be transient and would most 
likely be observed during their quiescent state at low lu- 
minosities. Muno et al. (2004) suggest that the faint hard 
source population (F < 1) can be consistent with mag- 
netic cataclysmic variable systems. It is suggested that 
the transient RLOF-fed NS-WD transients in quiescence 
studied in this paper may contribute to the faint popula- 
tion characterized by softer spectra. 

For the bright systems, 15 XRBs are found in our sim- 
ulated population of the MM03 field. These sources are 
characterized by luminosities ~ 10"^^ — 10"^^ ergs s~^ and 
are RLOF persistent and RLOF transient systems in out- 
burst. However, only a few bright transients have been 
observed in the field of the MM03 survey (M.Muno, pri- 
vate communication) , suggesting that the simulation leads 
to an overproduction of such systems (see §4) . 

For the WGL02 field, our standard model calculation 
shows ~ 60 X-ray sources brighter than 10^^ ergs s~^. 
Most of these systems in our simulation are persistent X- 
ray sources. However, only ^ 20 bright sources were ob- 
served (WGL02) with most as transient systems in out- 
burst in the Galactic Center. This suggests that the 
parameters intrinsic to our standard model affecting the 
bright population require modification. For the alternative 
models (full AIC kicks, lowered CE efficiency, reduced stel- 
lar winds and AIC not allowed) it is found that number of 
the bright sources is reduced to less than ~ 20. We note 
that the number of transient sources exceeds the number 
of persistent sources only for the simulation with reduced 
stellar winds. 

To account for the entire bright Galactic field XRB pop- 
ulation we increase the above number by factor of ~ 100, 
yielding ^ 6000 bright systems (4000 persistent and 2000 
transients in outburst). This is an overestimate, since only 
150 of these systems are observed. This overproduction of 
bright systems is common to all our models (including the 
neglect of AIC) and the other studies of the RLOF X-ray 
sources in the Galaxy (e.g., Pfahl, Rappaport, & Podsiad- 
lowski 2003) and implies an overproduction of the bright 
population in the Galactic center as well. We note that 
the overproduction may not be as great for the faint pop- 
ulation since the relative number of faint to bright sources 
can change from model to model. In particular the num- 
ber of bright sources decreases by a factor of 5, whereas 

^numbers from M. Muno, foreground sources excluded; private cc 



the number of faint sources is relatively unchanged for a 
model with reduced stellar winds. 

3.5. Comparison with other studies 

Willems & Kolb (2003) estimated the number of prc- 
LMXBs with main sequence companions less than about 
2 M0 in the Galaxy. Their number may be directly 
rescaled to the smaller expected population in the Galactic 
center, yielding ~ 10^ — 10^ sources for the WGL02 sur- 
vey field (B. Willems, private communication). Since the 
study of Willems & Kolb (2003) assumes no bolometric 
correction, we use our model with no bolometric correc- 
tion for comparison. In this case we expect 250 wind-fed 
low-mass X-ray binaries (prc-LMXBs), a number which is 
consistent with the Willems & Kolb (2003) findings. If 
the bolometric correction was applied, the expected num- 
ber of pre-LMXBs with luminosities ^ 10^° erg would 
decrease by a factor of ^ 2. 

A direct comparison with PRP02 is not easily accom- 
plished since the initial input (e.g., IMF, eccentricities) 
and treatment of the RLOF phases as well as model as- 
sumptions (e.g., natal kicks) difii'er. Despite the model 
differences, it is found that the numbers of wind-fed sys- 
tems predicted by PRP02 and in the present study for the 
WGL02 survey can be comparable (see below). In addition 
the X-ray luminosity distributions show that the majority 
of the synthetic wind- fed sources of both studies are char- 
acterized by rather low-luminosity ^ 10'^^ — 10^'^ ergs 
s~^ (although we are unclear about the adopted choice for 
the bolometric correction in PRP02). The primary differ- 
ence stems from the fact that the wind-fed sources in our 
study are characterized by low- and intermediate-mass MS 
companions, while PRP02 accounted only for the wind-fed 
sources with either intermediate- or high-mass compan- 
ions. 

In the present study the intermediate mass X-ray bina- 
ries (IMXBs) greatly dominate over HMXBs, in contrast 
to that found by PRP02. In fact, PRP02 obtain in their 
model K2 about equal numbers of HMXBs and IMXBs. 
This arises from their specific choice of kick scenario, which 
favors the survival of HMXBs. We prefer to use the ob- 
servationally derived distribution of kicks (Arzoumanian 
et al. 2002). With the use of rather standard kick veloc- 
ity distribution (model Kl), the PRP02 calculation results 
in the dominance of IMXBs over HMXBs (by a factor of 
about 5-10; as one can infer from the middle panel of their 
Figure 1). This result is already comparable to the results 
obtained in the present study. Moreover, there is an addi- 
tional factor which lowers the number of HMXBs as com- 
pared to IMXBs within our population synthesis model. 
Spec;ific;ally, PRP02 use a constant escape wind speed in- 
dependent of mass (or spectral type). However, it is known 
that this speed depends on a mass of star (e.g. Lamers, 
Snow & Lindholm 1995). For the mass range used in our 
study it varies by about an order of magnitude. The cap- 
ture rate (and thus resulting X-ray luminosity of wind fed 
system) depends strongly on this velocity. The higher the 
wind velocity, the smaller the capture rate and the lower 
the X-ray luminosity. In our simulation the most massive 
stars have wind speeds as high as 5 times the escape speed 
(Lamers et al. 1995), while PRP02 considers wind speeds 



7 



not greater than twice the escape speed. This effect results 
in a smaller number of bright X-ray binaries with massive 
companions (HMXBs) as compared to PRP02 study. 

For the WGL02 survey field PRP02 obtain a total of 
250 and 600 wind-fed neutron star systems with main se- 
quence companions more massive than 3 Mq for their kick 
models Kl and K2. The number of PRP02 systems is re- 
duced due to the detection limit appropriate for WGL02 
survey {L^ = 10'^^ ergs s~^) by a factor of 2.5 for a wind 
speed equal to the escape velocity or by a factor of 20 for a 
wind speed equal to twice the escape velocity. This leads 
to a reduction in the number of systems in the PRP02 Kl 
model from 250 to 12-100 systems in the WGL02 field. 

Since the Kl kick distribution of PRP02 (a single 
Maxwellian with a = 300 km s^) is comparable to the 
present study (bimodal distribution with 40% kicks with 
a = 90 km s~^ and 60% with a — 500 km s^^ - Arzouma- 
nian ct al. 2002), we compare our results with the results 
from model Kl. Our total number of wind- fed sources 
(standard model) for WGL02 field is 224 (see Table 1; 
note an increase by factor of 4 for the WGL02 field). In 
the following we give the different corrections factors in 
order to compare our results to that obtained by PRP02. 
First in the model with no bolometric correction (assumed 
to be the case in PRP02) we expect a greater number of 
wind- fed sources (320). Of these, only a 44 have luminosi- 
ties in excess of 10^'^ ergs s~^. Those systems with MS 
companions which are more massive than 3 Mq amounts 
to 32. Therefore, we conclude that the overall number of 
wind fed neutron star systems systems in PRP02 study 
for model Kl is in approximate agreement with the re- 
sults obtained in our study for similar model assumptions. 
However, PRP02 increase the number of WNS by using 
their alternative kick model (K2) , which can reproduce the 
entire population of X-ray sources in the WGL02 field. 

4. DISCUSSION 

The population of low luminosity X-ray binaries in the 
Galactic center region has been investigated within the 
framework of a population synthesis technique. Although 
it has been found that this population may be more het- 
erogeneous than previously predicted, our results suggest 
that accreting neutron star systems are not likely to be the 
major contributor to the faint X-ray source population in 
the Galactic center. 

The RLOF systems investigated in this paper provide 
an additional low luminosity component. Many of these 
additional systems are composed of neutron stars with 
white dwarf companions in ultra-short orbital periods 2 
hrs) in which mass is transferred at rates ranging from 
10^^^ — lO"""^" Mq yr^^ as a result of the action of angular 



momentum losses associated with gravitational radiation. 
Because the mass transfer rates are low, these systems are 
expected to exhibit transient quiescent/out burst behavior. 
For a small duty cycle of 1% it is expected that most of 
these sources would be found in their quiescent state with 
X-ray luminosities in the range of ^ 10^° — lO"^"^ erg sec~^. 

The population of hard X-ray sources (10^'^ — 10^^ erg 
sec~^) discovered in WGL02 survey was interpreted in the 
context of wind- fed NS-MS binary systems (PRP02). We 
find that such systems can be a contributor to the X-ray 
population in this luminosity range provided that rj^oi ^ 1. 
However, our results cannot explain the entire population 
since the predicted number of such systems (^ 30) is much 
lower than the number of observed sources (^ 100). 

The deeper MM03 exposure revealed a large population 
of faint {Lx ~ 10^° - 10^^ erg sec~^) point sources. The 
majority of these sources are characterized by hard spec- 
tra, and are too numerous to be explained by wind-fed 
NS/BH systems and too hard to be explained by RLOF- 
fed NS/BH systems. Muno et al. 2004 have suggested 
that the intermediate polar class systems can explain both 
their spectral properties and their observed number. The 
nature of the population characterized by softer spectra 
is unknown, although a number of candidates have been 
proposed. We have found that the number of RLOF-fed 
NS/BH transient systems in quiescence could be as high 
as 400 in comparison to the 652 observed, but this leads 
to an overproduction of bright LMXBs in the Galaxy. The 
sensitivity of our results to input parameters have been ex- 
plored, which can lead to a reduced contribution of RLOF 
sources and to an overall reduction in the number of syn- 
thetic soft NS/BH X-ray sources to a level further below 
that observed in MM03 survey. 

The identification of the faint soft X-ray sources will be 
essential for confirming their nature. Such studies may be 
fruitfully carried out in the low extinction regions near the 
Galactic center where the contribution from the low energy 
component of their spectra (corresponding to less than 1 
keV) may be detectable. At these lower energy ranges, the 
X-ray luminosities of our proposed candidates could be as 
high as - 10^^ - 10^3 ergs s-\ The Ghandra and HST 
surveys of Baade's window have been carried out (Grind- 
lay et al. 2003), and the forthcoming results may have 
bearing on the importance of this population. 

We would like to thank M.Muno, D.Wang, 
T.Maccarone, B.Willems, P.Podsiadlowski and N.Ivanova 
for useful discussions on the project and the anonymous 
referee for the comments which have helped to improve 
our study. This research was supported in part by the 
NSF under Grant No. AST-0200876 to RT and by the 
KBN Grant PBZ-KBN-054/P03/2001 to KB. 
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Fig. 1. — X-ray luminosity distribution for RLOF-fed (thick line) and wind-fed (thin line) binaries for the standard model calculation 



described in (§ 3.2). Note the very clear dominance of RLOF-fed systems in the faint luminosity regime {Lx 
simulated observed sample of Galactic center X-ray binaries. 
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Fig. 2. — X-ray luminosity distributions for RLOF-fcd (thick line) and wind- fed (thin line) binaries for different assumptions on the initial 
mass distribution of binary components (top panel) and common envelope efficiency as compared to our standard model. For details of 
alternative model assumptions see text (§ 3.3). Note the change in vertical scale on the panels. 
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Fig. 3. — X-ray luminosity distributions for RLOF-fed (thick line) and wind- fed (thin line) binaries for different assumptions concerning 
accretion induced collapse: full natal kicks incorporated (top panel) and AIC systems removed from the sample. Details on the models in 
§3.3. 
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Table 1 

X-RAY Binaries in the Galactic Center 



Type^ 


Wind-fed 


RLOF-fed 


Total: 

NS accrctor 
BH accrctor 


11.4% (56) 
11.0% (54) 
0.4% (2) 


88.6% (435) 
74.1% (364) 
14.4% (71) 


MS donor 
giant donor 
WD donor 


10.2% (50) 
1.2% (6) 
0% (0) 


21.6% (106) 
7.7% (38) 
59.3% (291) 


transient^ 

ultrashorf^ 
SN system^ 


0% (0) 
0% (0) 
11.2% (55) 


85.9% (422) 
60.7% (298) 
14.4% (71) 



^ X-ray active binaries (L^ > 10^^ erg sec~^) are listed. 
The relative occurrence frequency is given, followed by 
the occurrence frequency of specific types of accrctors and 
donors (MS main sequence star. WD white dwarf. NS- 
neutron star, BH- black hole). In parenthesis we list the 
expected number of sources in the Galactic center region 
corresponding to the Muno et al. (2003) survey (see § 3.4 
for details). 

''All the transients are listed here. For the assumed 
duty cycle of 1%, only 3-5 of these systems are expected 
to be observed in the outburst state. 

'^Systems with orbital periods shorter than 2 hrs. 



■^Systems in which there was a SN explosion. 



